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Abstract. This paper considers the segregation of a granular mixture in a rotating drum. Extending a recent kinematic 
model for grain transport on sandpile surfaces to the case of rotating drums, an analysis is presented for radial segre- 
gation in the rolling regime, where a thin layer is avalanching down while the rest of the material follows rigid body 
rotation. We argue that segregation is driven not just by differences in the angle of repose of the species, as has been 
assumed in earlier investigations, but also by differences in the size and surface properties of the grains. The cases of 
grains differing only in size (slightly or widely) and only in surface properties are considered, and the predictions are in 
qualitative agreement with observations. The model yields results inconsistent with the assumptions for more general 
cases, and we speculate on how this may be corrected. 

PACS. 47.55. Kf Multiphase and particle-laden flows - 83.70.Fn Granular Solids 

1 Introduction segregation or de-mixing, and even the parameters that control 

mixing and segregation are not fully understood. 
The production of many goods, ranging from pharmaceuticals 

and foods to polymers and semiconductors, depends on reli- Rotating drums, or kilns, are employed in industry to carry 

able uniform mixing of granular materials. Although there have out a range of oper ations; some examples are the calcination 

been several recent advances, particulate mixing is poorly un- of limes tone, reduction of oxide ore, clinkering of cementi- 

derstood, and one cannot a priori predict the effectiveness of tious materia i s> waste incineration and calcination of petroleum 

any mixing process. Indeed, mixing operations often result in coke 0wing to its industrial i mporta nce the rotary kiln has 
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for certain industrial processes, there is considerable evidence 
of segregation when the charge is a granular mixture of dif- 
ferent properties. A mixture of grains differing either in size 
and/or roughness [n| -\m, or density [ph when rotated in a drum 
is seen to undergo radial segregation, characterized by the for- 
mation of a core of smaller or rougher grains. This radially 
segregated core is seen to evolve into alternate bands of larger 
p ; __~ (or smoother) and smaller (or rougher) grains along the length 

Layer 

of the drum MM]. Axial segregation is found to occur over a 




Fig. 1. Streamlines in the active and passive layers in the rolling considerably larger time scale than radial segregation, which is 

regime of bed motion in a rotating drum. Grains cascade rapidly down usua u y complete within a few drum rotations. In a rotary kiln, 

the active layer, which is usually very thin in comparison with the radial se g re gation could lead to poor contact between the gas 

" ^ ' flowing above and the particles in the core, which would re- 
sult in poor heat transfer and/or lower rate of reaction. Axial 
provement in kiln performance may be achieved by better un- 

segregation, would lead to products of fluctuating qualities, 
derstanding grain transport during operation. 

In a rotating drum several regimes of bed motion have been Analysis of granular segregation has been undertaken in 
identified, namely, slipping, slumping, rolling, cascading, cataract-some recent studies via discrete element simulations [[§, and 
ing and centrifuging [f|. The most desirable bed motion for using coarse-grained continuum models ||-[[I(J. A recent ap- 
many industrial operations is the rolling mode, as it helps in proach, which is the one followed in this work, is based on the 
promoting good mixing of the particles along with rapid re- theoretical formalism of Bouchaud, Cates, Ravi Prakash and 
newal of the exposed material. The rolling regime is charac- Edwards [|TJ. Henceforth referred to as BCRE in this work, 
terized by two distinct regions: the "active" and "passive" lay- this study models grain transport on a sandpile surface, and 
ers [g]. Grains in the passive layer execute rigid body rotation has been employed in [^ for studying segregation during the 
along with the drum and the streamlines are circular (see fig- filling of a silo with a mixture of two species. The "minimal" 
ure 1). The active layer is usually very thin in comparison with model of ^ describes the case of grains with equal sizes but 
the extent of the passive layer. Grains roll down rapidly in the with small differences in surface properties. The same formal- 
active layer and streamlines have been observed to be straight ism was used in [^ to demonstrate the possibility of complete 
l]3| . segregation in a mixture of large smooth grains and small rough 

While rotating drums are used to mix particles of varying grains, and spontaneous stratification, i.e. alternating layers of 

sizes and shapes to obtain a homogeneous mixture necessary the two species, for a mixture of large rough grains and small 
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smooth grains. The minimal model of [|l2| was generalised by mobile grains in the passive layer at the interface 4> a (x, t) and 



[ 14 1 to accommodate grains differing slightly in surface prop- <j>p (x,t). The surface of the active layer is inclined at an angle 

erties as well as size. 6 r , the dynamic angle of repose), with the horizontal. 

In this work, we extend the formalism of Q to address the The conservation equations for the two species in the active 

problem of segregation in rotating drums of a granular mixture layer are 
whose constituents differ in size and surface properties. The 



problem of segregation in rotating drums has been addressed 



d t R a (x,t) = -d x (v(x,t)R a (x,t)) + T a 



in [|L5| for some special cases; however, our approach differs +d x (Dd x R a (x,t)) (1) 

from theirs in some fundamental ways, which we will point d t Rp{x,t) = ~d x {v(x,t)Rp(x,t)) + Tp 



out in this paper. 



+d x (Dd x R (x,t)) (2) 



o The model where v(x, t) and D are the mean velocity and the diffusivity of 

the grains, respectively. The active layer is assumed to be very 



Following BCRE, we assume the existence of a sharp interface 
between the active and passive layers. Grains in the active layer 
are referred to as rolling grains and those in the passive layer as 
immobile. Collisions between rolling and immobile grains oc- 
cur at the interface, which lead to exchange of grains between 
the two regions. We are interested in studying the case of a bi- 
disperse mixture consisting of a and (3 type grains, which differ 
in surface properties and/or size. 

We call R a (x,t) and Rp (x, t) the number of a and /3 type 
rolling grains in the active layer per unit length of the inter- 



thin, and variations across it are not resolved in this model. All 
the interesting physics in this model is buried in r a and Fp, 
which are the rate of conversion of immobile grains of each 
type into rolling grains by collisions at the interface, and we 
will specify their functional form shortly. 

In the passive layer, grain transport is only by advection due 
to the rotation of the drum. The equation for conservation of a 
type grains, integrated across the passive layer, is 



Pnd t / 4> p a {x,y,t)dy = cf> a p n v g .n- T a , (3) 

face, and R(x,t) = R a (x,t) + Rp(x,t) the total number of 

rolling grains (per unit length) in the active layer. The local and similarly for (5. Here, ^,(x,y,t) is the number fraction of 

slope of the interface is 8(x,t), the thickness of the active layer a type grains inside the passive layer, vg is the velocity field 

is S(x,t) and the height of the passive layer is h(x,t) (see fig- due to solid body rotation, p n is the number density of grains 

ure 2). As the material in the passive layer only undergoes solid in the passive layer (taken to be constant), and n is the unit 

body rotation, it is sufficient to specify the the concentration at normal to the interface, pointing from the passive layer to the 

the top surface to fully specify the concentration throughout active. (Due to the circular streamlines in the passive layer 0^ 

this layer. Consequently, we define the number fractions of im- can be related to the number fraction at the interface <fi a .) The 
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Fig. 2. The coordinate frame for our model equations, showing some 
of the variables of interest 

conservation for all grains then assumes the form 



p n d t h = /9„v .n- ^ r i 



(4) 



2.1 Model for r, 

As in BCRE, we write an expression for F a in terms of the 
rolling grain densities Ri and the number fraction of immobile 
grains </>, at the interface. Boutreux and de Gennes (1996) con- 
sider the four possible outcomes of binary collisions between 
rolling and immobile grains at the interface; amplification, cap- 
ture, exchange and recoil. Amplification of species a refers to 
the event of an a type immobile grain being dislodged into 
the active layer by a rolling grain. This could happen either 
by auto -amplification, when the colliding grain is also of type 
a, or by cross-amplification, when it is of type j3. The auto- 
amplification (cross-amplification) rate of a type species will 
increase with R a (Rp), 4> a and the collision frequency /. Cap- 
ture of species a refers to the event of an a type rolling grain 



being rendered immobile and absorbed by the passive layer. 
Arguing similarly for the rate of capture (and noting that ex- 
change and recoil leave the net transfer of grains unchanged), 
we write the simplest form for r a as 



r a =i a aa ^ r R a <i> a g a ^) + i a a p^ r R^ a g a W 



-j c a ^R a g c W 

dp 



(5) 



wherein we have taken the collision frequency, /, to vary as 
-T-, dp being the average grain diameter. The coefficients j^ a , 
7^ o and 7^ are constants of order unity and g a (ip) and g c (i/j) 
functions of the slope of the interface. 

We note that unlike in previous studies investigating granu- 
lar segregation Q12, 14 15|], we do not set the velocity v to a con- 
stant, but allow it to vary along the length of the active layer by 
solving the momentum balance in conjunction with the species 
balances ([!]). 

To complete the formulation, the functional forms of g a (i/j) 
and g c (ip) must be specified. Recognizing that the rates of am- 
plification and capture depend not just on the local slope of the 
interface but on the difference between the local slope and the 
angle of repose, we set ip = 9(x) — 9 r (x), where 9{x) is local 
slope of the interface and 6 r (x) is the angle of repose for the 
mixture. Next, J 1 ™ (P°) should rise (diminish) with ip, but we 
do not expect it to vanish altogether at i\) = 0, since we expect 
a fraction of the grains striking a surface at the angle of repose 
to dislodge other grains (get captured). However, for large neg- 
ative (positive) values of ip, the rate of amplification (capture) 
should become negligibly small. In accordance with these ar- 



Saikat Chakraborty et al.: Analysis of radial segregation of granular mixtures in a rotating drum 5 

guments, we assume the following continuous forms: through the slope functions (eq. ||). Our results, shown in sec- 



. , . tions 3.1 - 3.3, illustrate quite clearly that the qualitative na- 

i -r- exp^yj) ture f s j ze an( j roughness segregation is different, which is in 

A 
9 ^' ~ i -)_ cxp(Btb) agreement with previous experimental observations P] - p]. 

d t We allow the the angle of repose to depend on the grain con- 

where A and B are constants of O(-r-) and O(-j-), respectively. 

centration in our model, but assume it to be constant in this 
paper for the sake of simplicity. 



The notable feature of the above forms is that g a + g c is con 



2.2 Momentum Balance 



stant - while we are unable to ascribe a physical meaning to 
this constraint now, it is essential if the solutions are to be sym- 
metric about x = 0. The issue of symmetry is commented on in 

greater detail in section EOl. Previous studies that used the same theoretical framework made 
We note here that MakselQ considered special cases of the unrealistic assumption that the velocity of grains in the ac- 
the above slope functions for very small and very large dif- tive is constant - clearly, the velocity must start with zero at the 
ferences of angles of repose. He also assumed that differences t0 P end of the active la y er ' reach a maximum at an intermed- 
in size and surface properties between the species result only in ate P oint ' and a 8 ain vanish at the lower end of the active la y er - 
differences in their angles of repose, which is the ultimate cause We allow the velocit y t0 var y alon 8 the len 8 th of the la y er ' and 
of segregation, and that the coefficients for auto-amplification write a momentum balance for the rolling grains, averaged over 
and cross-amplification ( 7 ° Q and 7^) are equal. We believe the thickness of the active layer: 

that ascribing segregation purely to a difference in the angle of Q t (# w ) + g x (jj^ =v V^ p. + Rg s [ n g r 

repose is unrealistic; it surely does not address segregation of 



grains of different sizes, but made of the same material. In this 
case, segregation arises because the probability of a larger grain 
dislodging a smaller one is greater than that of the larger one 
dislodging another larger one, as we shall demonstrate below. 
We propose that segregation due to size and surface prop- 
erties are driven primarily by amplification and capture respec- 



—<Tf - cr c , (8) 

The first term in the right hand side accounts for x momen- 
tum influx due to conversion of immobile grains into rolling 
grains, the second is the momentum generated due to the gravi- 
tational body force, and the third and fourth terms are the shear 
stresses at the interface due to frictional grain interactions at 
the interface and grain collisions [n6h, respectively. We assume 
a Coulombic form for the frictional stress, 



tively, as elaborated in section 3.1. This has been done by al- 
lowing (unlike [|l5||) the coefficients 7^ Q , 7%, 7° etc. to de- 
pend on the size, roughness and type of interacting grains. In 
addition, the effect of differences of angle of repose caused by where fj, is the coefficient of friction of the material. The col 
differences on size or surface properties has been incorporated lisional shear stress, a c , is generated by collisions of the ex 



Of = uRg cos 9 r , 
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change and recoil type, discussed earlier, that do not lead to with, 
inter-conversion of grains but do result in loss of momentum. 
Therefore, we choose a constitutive form of a c that is similar 
to the momentum flux due to inter-conversion, 



r a = -f aa vR a (j) a f a (^) + 7 Q/3 v R/3 0q f a ( ip) 

- 7 >i? Q / c (V) (10) 



<r. = A*£, rW = ^W ;rW = — l (11) 

d p 1 + cxp(i/)) 1 + exp(V') 



where A, like A, is a constant of proportionality of 0(-j£). where 



2.3 Scaling and Leading Order Analysis 

To clearly illustrate the relative magnitude of the different terms 
in our model equations, we define the following dimensionless 
variables, 



- ^ 
dx 



The above forms of the slope functions result upon setting AL/d p 
Bdp/L = 1, which we have done for simplicity. 

We then consider the case of small e, as is the case in most 
practical situations, and seek the to leading order solution. In 



x i) o R h 

./■ = — , u* = , S* = — , Ra — — — , h* = — this approximation, the governing equations reduce to the fol- 



lowing set: 



L uj L dp p n L L 

and the following dimensionless parameters, 

d(vR) v- r nr> 

e ~ r > ' "~ r ' "^ "~ . r j ' « — 2 r ' A ~~ j i=<*,0 



L L uiLdp uj L d n 



d{vRa) = r Q , (13) 



where L is the half-length of the free-surface, lu is the rotational dx 

velocity of the drum, d p is the average grain size, r is the radius 

- x = Yl r *> < 14) 

of the drum and p„ is the number of grains per unit area of the i=a,/3 



passive layer. 

The active layer is expected to be only a few grains deep, 



i X 1 (x , H-^/ 



and its thickness 6 is therefore scaled with d p . The diffusion ,,_ 2 \ 

— = v 2_. ^i + R,9 [sin^r — M cos9 r ] 

coefficient, D, being a material property scales as the product "-x i- a a 



of the velocity and the grain size. Using the above normalized 
variables, and dropping the asterisks for the sake of simplic- -^ 
ity, (|l]) tranforms to the following dimensionless form at steady 
state: 



XRv 2 , (16) 



R a + Rp = R, (17) 



d(vR a ) d 2 R 

K dx ' = eD-^-j + T a (9) < f >a + ( j )ft = i. (18) 



Saikat Chakraborty et al.: Analysis of radial segregation of granular mixtures in a rotating drum 

An additional condition is required for solving for ip a , and Substituting ( p0[ ) and ( p"5[ ) in ( pj| ) and simplifying yields 

this is the loading condition specifying the total amount of ma- d u (u — 6 )x 



dx R v 

terial loaded into the drum, 



(22) 



L j-y a where u a = R a /R is the number fraction of a grains in the 

/ <%dydx = V$ a , 
L-'Vd active layer. Now (|14|), and ( |15| ) can be used to express <fi a (x) 



where V and <P a are the total volume and number fraction of a in terms of Wq and obtain the solution of Q The constant of 



type grains in the entire charge, respectively. Since ^ a is con- integration for this solution must be evaluated by enforcing the 

stant along each circular streamline in the passive layer, the i oac ji n g condition given by dl9|) 

above condition reduces at leading order to the following di- We consider various cases where the grains of the two species 

mensionless form. differ either in size or in surface properties. The mean number 

r 1 

/ 6 (x)xdx = <P (19) fraction <f> a has been set to 0.5 in all the cases. 



3 Model Predictions 



3.1 Grains differing only in surface properties 



_ , . . r—i . .rrri , , . , , , ... It is reasonable to suppose that amplification is largely guided 

Substituting fll4| ) in (|12J), and applying the boundary condition 

by differences in sizes while capture is driven by differences in 



vR = at x = — 1, we get 



vR=£^l. (20) 

2 



surface properties. Then for this case, we may set 

~; a ~, a ~A ~A ~A 

laa = lafi = Ipa =1 013=1 ■ 



Substituting (EOb in (161), and solving for the velocity subject to 

the boundary condition « = at i = -1 we get The differences in surface properties would, of course, lead to 



where 



(j the capture coefficients being different, and we set 

1% = 7 C , l a = si c p . 



-(l-exp[-A(l + x)]), (21) 



C = g[sa\9 r — ficos9 r ]. 

The above velocity profile does not vanish at the lower end 
of the active layer, y = 1; this boundary condition may be en- 
forced by putting in a bulk viscosity for the flowing medium. 



A value of s > 1 implies that a type grains are rougher than 
/3, and hence more easily captured by the passive layer. With 
these simplifications, it follows that the solution of (|13|)-(|l5|) 

is: 



U a „ ,, „2\s-l 



and — =C 1 (l-x 2 Y 



The correction to the velocity would be of higher order in e Va ~ \ _|_ ( s _ l)u a (1 — u a Y 

everywhere, except in a thin boundary layer near y = 1. Also (23) 

notable is that while the velocity is not symmetric about x = 0, where C\, the constant of integration is obtained by numerical 

the flux v R : is. integration of (|l9|). Figures 3 and 4 present the variation of the 
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Fig. 3. Variation of the number fraction of a type grains (the rougher 
species) in the active layer along the length of the interface for the case 
of grains differing only in surface properties. The relative roughness 
of the a species increases with the parameter s. 

number fraction of a type grains in the active and the passive 
layers, respectively, along the length of the interface for two 
cases, one where the a type grains are slightly rougher than the 
(3 type grains (s = 1.5), and another where the a type grains 
are much rougher (s = 10). The model predicts depletion of 
the rougher grains at the two ends of the drum and accumula- 
tion near the center of the drum. The degree of segregation, of 
course, depends upon the difference between the surface prop- 
erties of the two species, quantified in this model by s. 

3.2 Grains differing slightly in size 



• -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 

X 



Fig. 4. Variation of the number fraction of a type grains in the passive 
layer along the length of the interface, for the same parameter set as 
in Fig. 3. 

course, that of constant number fraction in the active and pas- 
sive layers, 



U a (x) = (f> a (x) = <P a 



(24) 



where <P a is the mean number fraction of a type grains during 
filling of the drum. For a mixture comprising species that differ 
very slightly in size, we allow the amplification constants for 
the two species to differ by a small amount -d. The probabil- 
ity of a large rolling grain dislodging a large immobile grain is 
equal to that of a small rolling grain dislodging a small immo- 
bile grain. Hence, 






(25) 



. Following our arguments in § 3. 1 , we will assume that the dif- 
ference in size leads to a difference in the rate of amplification 

of the two species but not that of capture. We investigate the However, if the (3 type grains are larger, the probability of a (3 

case where grains have identical surface properties but differ type grain dislodging an a type grain is greater than that of it 

slightly in size. For the case of grains differing neither in size dislodging one of its kind, and the probability that an a type 

nor surface properties (i.e. identical species), the solution is, of grain dislodges a /3 type grain is less than that of it dislodging 
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one of its kind. Hence, for small differences in sizes, we set 



7^=7 Q (1 + ^), 



7^=7 a (l-tf)- 



(26) 



(27) 



We now seek the solution as a perturbation to the uniform so- 
lution (g4j) to leading order in #, 

u a (x) =$ a + #u' a (x), 4> a {x) =$ a + 0<fJ a (x) (28) 



and solve for u' a (x) and <f>' a (x). Substituting ( |25| ) - < j28j ) in (|1 
- (O), we can express (j)' a (x) in terms of u' a (x) and x as 

fax) = u' a (x) - 7 C ^^( £* V ~ X a )|~- (29) 

7 Rv + 7 x Rv 

Using (£8|), we can write ( p2[ ) as: 

du' a = « ~<j>' a )x 
dx Rv 

We now consider special cases for which simple analytical 

solutions of (f30h are possible. The first is the case of the grains 

of both species being very smooth, i.e -Q -C 1. The solution is 




Fig. 5. Variation of the number fraction of a (smaller) and f3 (larger) 
type grains in the active layer along the interface. For this case, both 
species are very smooth, and the parameter $ is set to 0.1. 



7 



u a (x) =$ a +d-^$ a $ r3 \n{l-x 2 ), 

7 

7 C 



(31) 



b a {x)=$ a + {>^$ a $p[\ + Hl-x 2 )]. (32) 

7 



Figures 5 and 6 present these results for # = 0.1, where it is 
clear that the central core is richer in smaller grains, in agree- 
ment with experimental observations. 

The second case is when both species are very rough, i.e 




Larger Grain 
- Smaller Grain 




0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 



-% 3> 1. The solution is 



7 



u a {x) =$ a -$$ a $ ln(l-x 2 ), 
M*) = *o - ®$ a $p [1 + Ml - x 2 )]. 



(33) 



(34) 



Figures 7 and 8 present the profiles of u a and <fi a for $ = 0.1, 
where it is apparent that, unlike in the case of smooth grains, 



Fig. 6. Number fraction of a (smaller) and j3 (larger) type grains in 
the passive layer as a function of position on the interface. Parameter 
values the same as in Fig. 5. 

the smaller grains are depleted near center of the drum and tend 
to accumulate at the two ends. This is due to the fact that if 
both species are equally rough, the larger grains will roll further 
down the active layer as their inertia is greater. No experimental 
results are however available to validate this prediction. 



10 
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Fig. 7. Profiles of number fractions of a (smaller) and j3 (larger) type 
grains in the active layer when both species are very rough, for ■& — 0. 1 
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Fig. 8. Profiles of number fraction of a (smaller) and /3 (larger) type 
grains in the passive layer at the interface, for the parameter set as in 
Fig. 7. 

3.3 Grains differing widely in size 

For large differences in sizes between the two species, experi- 
mental observations in sandpiles reveal either complete segre- 
gation or stratification; it has been reported [|l7|| that this occurs 
when the size ratio exceeds 1.5. Boutreux et al. [|l4|] proposed 



that size segregation is then due to percolation or "kinematic 
sieving", as the smaller rolling grains tend to fall through the 
gaps between the large grains, forming a sub-layer above the 
interface. The larger rolling grains are therefore not in contact 
with the bulk, and are not captured to the extent they would 
have been if in contact. To account for this phenomenon, we 
follow [ |l4[ | and modify r a by a factor exp(pu a ), which mim- 
ics the screening of the large grains from the interface by the 
sub-layer of the small grains. The dimensionless parameter p, 
which is proportional to the size ratio of (3 to a type grains, 
measures the degree of percolation. 

For a large difference in sizes, segregation due to percola- 
tion effect would far exceed the segregation effects due to dif- 
ferences in amplification and capture coefficients, as suggested 
by [|l4|]. Hence we assume, 



7 Q q = lap = Ifia =l0p=l 



(35) 



and 



7 Q = 7/3 = 7 ' 



(36) 



Using ([35]) and @, we solve Q - @ to obtain: 



<l>a(x) 



u a (x) 



u a (x) + (1 - u a (x)) exp[pu a (x)] ' 



(37) 



and solve ( |22| ) numerically. The results shown in figures 9 and 
10 for p = 2 indicate that the larger grains are heavily de- 
pleted in the central region of the drum and accumulate at the 
ends both in the active and passive layers. This is in qualitative 
agreement with experimental observations am. 



Saikat Chakraborty et al.: Analysis of radial segregation of granular mixtures in a rotating drum 



11 




M ^ ¥ is X 2S ^ X ?£ 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.E 

a; 



Fig. 9. Profiles of the number fraction of a (smaller) and (3 (larger) 
type grains in the active layer for the case of grains differing widely in 
size. See text for parameter values 
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Fig. 10. Profiles of the number fraction of a (smaller) and f3 (larger) 
type grains in the passive layer at the interface for the case as in Fig. 
9. 

3.4 Other cases: Model inconsistency 

The solutions described in the previous section are symmetric 
about the mid-plane x = 0. However this is not a generic re- 
sult, rather only for the special cases we have considered, and 
also because the slope functions have been chosen to satisfy 
g a + g c — constant. If the grains differ in size and shape, or 



if the slope functions do not satisfy the above constraint, the 
solutions for Ri and fa are not symmetric. This result is incon- 
sistent with our assumption that the passive layer executes solid 
body rotation: circular streamlines in the passive layer necessi- 
tates symmetry of Ri and a>i 

It thus appears that the assumption of solid body rotation 
in the passive layer must be relaxed for consistency, and so- 
lutions in general need not be symmetric about the mid-plane. 
This may be accomplished by postulating an interface of finite 
thickness within which streamlines are not circular, which is 
the approach we are following currently and hope to report the 
results soon. The recent study of Shinbrot et al. [ |l8| ] reports the 
existence of an interface between the active and the passive lay- 
ers, wherein deformation occurs by periodic stick-slip motion. 
In this light, further experimentation is necessary to ascertain 
whether the concentration of species in the active and passive 
layers are symmetric. 

4 Conclusions 

Granular mixtures, consisting of grains that differ in size and/or 
surface properties, have been observed to undergo segregation 
when they are either poured in a heap or rotated in a drum. 
Several analytical models have been proposed to predict segre- 
gation in case of a two dimensional sand-pile. The theoretical 
formalism proposed by Bouchaud et al. [[FT]] and later extended 
by pi2[pj| ,p4j], successfully addresses the problem of segrega- 
tion in poured heaps. The present work attempts to extend this 
framework to the rotating drum problem. 

We argue that the qualitative nature of segregation of two 
mixtures, one of which differ in surface properties and the other 
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provides a canonical framework under which segregation of 
grains differing in size and shape could be studied in various 
geometries. The model gives satisfactory predictions for grains 
differing either in surface properties or size. These results are 
in good qualitative agreement with experimental observations. 
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it. 
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